Magnetic control of the interaction in ultracold K-Rb mixtures 
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We predict the presence of several magnetic Feshbach resonances in selected Zeeman sublevels of 
the isotopic pairs 40 K- 8r Rb and 41 K- 87 Rb at magnetic fields up to 10 3 G (1G = 1CT 4 T). Positions 
and widths are determined combining a new measurement of the 40 K- 87 Rb inelastic cross section 
with recent experimental results on both isotopes. The possibility of driving a K-Rb mixtures from 
the weak to the strong interacting regime tuning the applied field should allow to achieve the optimal 
conditions for boson-induced Cooper pairing in a multi component 40 K- 87 Rb atomic gas and for the 
production of ultracold polar molecules. 
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The possibility of controlling the interactions in an ul- 
tracold atomic gas using magnetically tunable Feshbach 
resonances Q] has been successfully exploited in recent 
remarkable experiments. It has indeed allowed to reach 
Bose-Einstein condensation (BEC) for species difficult to 
condense in zero field and to study their stability prop- 
erties 0, Q , to form bright solitons Q , to produce ultra- 
cold homonuclear molecules [j| and a strongly interact- 
ing Fermi gas . The presence of Feshbach resonances 
between atoms of different species could open new pos- 
sibilities, thereby in this paper we carry out a detailed 
investigation in this direction. 

We focus on two potassium-rubidium mixtures that 
have allowed the production of a two-species BEC ( 41 K- 
87 Rb) 0,0] and of a strongly interacting Fermi-Bose mix- 
ture ( 40 K- 87 Rb) 0. We predict the positions and the 
widths of several Feshbach resonances in selected mag- 
netic sublevels of the ground electronic state by combin- 
ing collisional measurements on the two isotopes. These 
resonances might for instance allow the observation of 
boson-induced Cooper pairing or the formation of 
ultracold bosonic or fermionic polar molecules 

We use for our calculations the collision model pre- 
sented in Ref. |l2| . where the Hamiltonian now in- 
cludes the Zeeman interaction with the external magnetic 
field 0| . Denoting "a" and "b" potassium and rubid- 
ium, collisions are labeled as (/ m/ )+(/iiW/i)), where 
f a , fb are the hyperfine angular momenta with projec- 
tion mf a , rrifi, along the quantization axis. Please recall 
that the model is parameterized by the singlet and triplet 
s-wave scattering lengths a s , t and by the long-range dis- 
persion coefficients C n with n=6,8,10. We adopt a ±2% 
uncertainty in the highly accurate van der Waals coeffi- 
cient C 6 = 4274 a.u. (1 a.u. = 0.0957345 xl(T 24 J-nm 6 ) in 
order to account for the combined uncertainty in all the 
C„'s. In the absence of weak relativistic spin-interactions 
the projection mf = nif a + nifb and the rotational an- 
gular momentum of the atoms about the center of mass 
£ are good quantum numbers for the collision event. In 



this work we limit ourselves to t = collisions which, 
at least in the absence of I > low-energy resonances, 
give the dominant contribution to the interaction in the 
ultracold regime. 

Two key measurements on 40 K- 87 Rb mixtures pro- 
vide the best range of input parameters to our model 
available to date. First, the elastic cross-section for 
(9/2, 9/2)+(2, 2) has been directly measured by inter- 
species thermalization |9j and by damping of coupled K- 
Rb dipole oscillations |14| . Combining the systematic 
uncertainties in these two determinations we estimate 
a t = a§§.22 = -410±|ga (1 a = 0.0529177 nm). Mass 
scaling of this value to 41 K- 87 Rb allows us to strongly 
tighten the upper bound on the triplet scattering length 
of Ref. [13 , leading to a* (41 -87) = I56±l° and to reduce 
the original uncertainty in the number of triplet bound 
states to iV b '(41 - 87) = 31 ± 1. 

Second, in order to determine a s we measure the 
(9/2, 9/2)+ (2, 1) atom-loss rate in a nondegenerate sam- 
ple. We prepare a doubly spin-polarized sample in the 
magnetic trap containing about 10 4 K atoms and 10 5 
Rb atoms at temperatures in the range (300 — 900) nK, 
and then transfer about 30% of the Rb sample in the 
(2, 1) state using a fast radio-frequency sweep. After 
a short rethermalization period we measure the time- 
evolution of the number of K atoms, which decays 
through inelastic collisions with Rb according to the 
rate equation n K (t) = -G inel nf /2 9/2 (i) n^(t). We find 
Qinei _ 5.0(2.5) xlO -12 cm 3 s -1 , where the uncertainty 
is dominated by that on the atom numbers. In this work 
we make a conservative choice and adopt only the order 
of magnitude G = (10~ 12 - lO" 11 ) cm 3 s _1 . 

Figure shows G mel numerically calculated as a func- 
tion of a s varying a t within the bounds determined 
above. The horizontal axis is labeled by a scaled "scat- 
tering phase" £ s = (2/7r) arctan(a s /a sc ), where a sc — 
^(2^iCe) 1 ^ 4 » 72ao is a typical length for a purely van der 
Waals potential, with /i the K-Rb reduced mass. To be 



cible even in the ultracold regime, see Ref. [if 




FIG. 1: The numerically calculated inelastic rate G me 
for 40 K(9/2,9/2)+ 87 Rb(2, 1) collisions at temperature T = 
500 nK as a function of the singlet "scattering phase" (see 
Text) for at = -330 ao (full line) and at = -490 ao (dashed 
line). The horizontal lines correspond to the bounds deter- 
mined experimentally. 



TABLE I: Compendium of K- Rb scattering properties. 
The zero-field s-wave elastic scattering length a is shown for 
different hyperfine sublevels. The uncertainty in the reso- 
nance positions -Bo results from the uncertainty in the model 
parameters quoted in the text. The resonance widths A are 
shown for the nominal values a 3 — — 185ao, a t = — 410ao and 
C(, — 4274 a.u. of the scattering lengths and of the van der 
Waals coefficient. 
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noted the strong suppression of the inelastic rate for a s ~ 
a t . Comparison with the experimental value leads either 
to a very large (positive or negative) or to a negative sin- 
glet scattering length a s (40 - 87) = -185+|25 a o- Mass- 
scaling with iV^(41— 87) = 98±2 gives a singlet scattering 
length for the bosons pair a., (41 — 87) = S.O+^ao and, by 
compatibility with the a s (41 — 87) of Ref. [12|, rules out 
the range of very large |a s |(40 — 87). With these a s , a t we 
can predict the zero-field scattering length and Feshbach 
resonances of different magnetic polarization states. Let 
us first consider the boson-fermion pair 40 K- 87 Rb. 

The inverted hyperfine structure of the fermionic 
potassium, with the / = 9/2 lying below the / = 7/2 
level, leads to a peculiar and very favorable feature of 
this pair, namely the stability under spin-exchange col- 
lisions of any combination where either atom is in its 
lowest Zeeman sublevel. We therefore find it important 
to summarize in Tab. [I] the scattering properites of the 
(9/2, — 9/2)+(l, nift) states and of at least the most eas- 
ily experimentally realizable (9/2, m/ a ) + (l, 1) magnetic 
states. Note all of these hyperfine states have a relatively 
large ( i.e. \a\ > a sc ) and negative zero-field scattering 
length. This ensures that the two clouds will remain mis- 



A magnetic Feshbach resonance arises when a molecu- 
lar bound-state crosses an atomic threshold [l^. Across 
a Feshbach resonance the scattering length a presents a 
typical dispersive dependence on the strength of the ap- 
plied field B, a(B) = a(l — B A B(j ), where B is the field 
strength on resonance, A is the resonance width and a is 
the off-resonance background scattering length. The res- 
onance parameters Bq and A have been extracted from 
an accurate close-coupled calculation on a fine i?-grid and 
are shown in Tab. [fl 

Homonuclear collisions between alkali atoms in a mag- 
netic field have been studied in great detail. Only a rel- 
atively small amount of Feshbach spectroscopy is often 
necessary as input to theoretical models in order to gain 
an impressive power in determining the interaction po- 
tential and in p redicting the position of additional reso- 
nances 0, Il8j . The parameters a s ,a t and C„'s of our 
collision model can in principle be determined following 
a similar procedure. At least three resonances between 
the ones shown in Tab. [I] would be needed as input to 
the model. Inclusion of the high-field (9/2, 7/2)+(l, 1) 
resonance would improve the procedure since this latter 



resonance is nearly purely triplet in character with av- 
erage electron spin (S 2 ) «1.95 and very sensitive to Cq, 
SB /B Q = -1.2 5C 6 /C 6 . 

The relatively strong and tunable interaction be- 
tween different magnetic polarization states opens up the 
possibility of realizing Bose-Fermi mixtures where the 
fermions undergo a Cooper-pairing phase transition to 
a superfluid state via exch ang e of density fluctuations of 
the background BEC 0, llflj - Due to identical particle 
symmetrization, direct as well as boson-mediated s-wave 
collisions between a pair of fermionic K atoms are allowed 
only in a Fermi gas with at least two internal components, 
which we denote as "1" and "2" in the following. We also 
let dRb be the intra-species scattering length for collisions 
of Rb atoms, cik the intra-species scattering length for 
collisions of K in state "1" and K in state "2" and a,KRb,i, 
with i = 1,2, the inter-species scattering lengths for col- 
lisions of Rb and K in state "i" . 

The critical temperature T c for boson induced s-wave 
Cooper pairing can be given an analytic expression valid 
in the limit of low boson-densities and of a BEC larger 
than the fermionic cloud, see Eq. (13) of Ref. 0. As 
long as a,KRb,ia<KRb,2 >> dRbdK, the critical tempera- 
ture is independent of cik, and increases exponentially 
with aRb/{kFaKRb,2<iKRb,i), where 1zf is the Fermi vec- 
tor. It must be noted however that as soon as kFa 2 KRb i ~ 
dRb, the z-th component of K and Rb will become dy- 
namically unstable and collapse [Til l20j| . As a con- 
sequence the optimal condition for Cooper pairing are 
aKRb,i ~ a>KRb,2 = 0-KR.b and a 2 KRb kp ~ Or&. 

From Tab. [i] a mixture composed of K(9/2, — 9/2), 
K(9/2, —7/2) and Rb(l, 1) seems to be the most appro- 
priate in order to fulfill this condition. However, for 
the nominal values of the zero-field scattering lengths 
between these components, we estimate that at least a 
factor of four increase in the average trapping frequency 
would be necessary with respect to our current exper- 
imental setup. If the actual axRb,% turned out to be 
smaller than the nominal value, this factor would increase 
even further. Fortunately, at least for some potential pa- 
rameters, a pair of resonances exist that allow simulta- 
neous tuning of axRb,i and a^Rba to a large value at 
the same time, making such a large trapping frequency 
unnecessary. 

A very favorable situation is shown in Fig. [2 where 
for B pa 725.89 G we get axRb ~ — 687ao . In order to 
have this value stabilized to within 10 ao the magnetic 
field must be stabilized to within 0.1 G, which is a strin- 
gent yet possible to be met requirement . The field at 
which aKRb,i — axRb,2 and the common value axRb are 
extremely sensitive to the details of the potential, and a 
conclusive word on what is the largest achievable com- 
mon aKRb is at this stage left to forthcoming precision 
experiments on K-Rb Feshbach resonances. As a trend, 
|et Sit | smaller than the nominal values tend to have a good 
working field for inducing the Cooper pairing instability, 
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FIG. 2: Field dependence of the 40 K(9/2, -7/2)+ 87 Rb(l, 1) 
(full line) and 40 K(9/2, -9/2)+ 87 Rb(l, 1) (dashed line) scat- 
tering lengths for a s — — 105ao and a t = — 380ao and the 
nominal Ca. The arrow shows the best working point for 
inducing the Cooper pairing instability (see Text). 

for \a s ,t\ larger than the nominal values (which however 
implies larger zero-field a's) such a field may not exist. 

We now consider the boson-boson pair 41 K- 87 Rb. We 
summarize in Tab. [[Tithe collision properties of the states 
(2, 2) + (2, 2) and (1, 1) + (1, 1), which are collisionally sta- 
ble for any field, and of the state (1,— 1)+(1,— 1) which 
is only stable up to a threshold value B l ss 143.45 G, 
where the (2,— 2) + (l,0) channel opens up and becomes 
the stable one. The 011,11 and ai_i.i_i, which are equal 
at zero field, are relatively large. 

The (1, 1)+(1, 1) collision always presents at least two 
broad and one narrow resonance below 10 3 G, see Tab.lTTl 
and Fig. [21 The analytical expression for the scattering 
length of a pair of overlapping resonances with positions 
Bq ± and widths A± 

<g), (g-go,+)A_ + (B-B ,_)A + 
a (B - Bq,+)(B - Bo,-) 

has been used to fit both low-field resonances in the ta- 
ble. The lowest-field resonance exists only if 011,11 > 0, 
for an, 11 < it turns into a real bound-state below the 
(1, 1) + (1, 1) atomic threshold and cannot be observed in 
a standard scattering experiment. 

The (1,— 1)+(1,— 1) collision presents a feature near 



TABLE II: Same than Tab.Q]but for 41 K- 8r Rb. The resonance 
widths A are shown for the nominal values a 3 = 5.0ao, at = 
156ao and Ce = 4274 a.u. of the scattering lengths and of 
the van der Waals coefficient. In cases where the sign of a 
is undetermined we have shown the nominal value, the upper 
bound for a < and the lower bound for a > 0. 
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the threshold value B l . In this situation virtual excita- 
tions to the (2, — 2)+(l, 0) channel can strongly modify 
the real part of the scattering length (see Fig.^J), defined 
as 7£ai-i,i-i = — 1Z lim^o Si-i,i— i/(2ik), with S the 
scattering matrix and k the wave vector of the incom- 
ing atoms. Such fluctuations- induced feature is always 
present in our range of parameters if > but 

turns out to be important for ax_i,i_i < only for very 

large |<xi 1,1—1 1 ■ To be noted the cusp at B = B t that 

follows from the analytical properties of the scattering 
matrix near the opening of a new threshold |2l| . 

Finally, a field induced resonance in the (2, 2)+(2, 1) 
inelastic loss-rate has been included in Tab. [H] since it 
may be important in order to constrain the collision 
model from Feshbach-resonance spectroscopy carried out 
on the boson-boson isotope. The width A shown in the 
table is in this case the inelastic one, defined through 
the expression of the total scattering eigenphase across 



resonance, 5(B) = do — arctan ■ 



with 60 the back- 



(B-Bo) ' 

ground scattering phase. Moreover, near this resonance 
we find a dramatic reduction of the (2, 2) +(2, 1) inelastic 
decay rate down to G » 10~ 15 cm 3 s _1 , a value that en- 
sures the collisional stability of a two-species (2, 2)+(2, 1) 
BEC. 

An interesting application of our results is to the pro- 
duction of polar molecules via photoassociation, as pro- 
posed in Ref. [llj . First step in this scheme is a Mott- 
insulator transition of a double K-Rb BEC into a state 
where only pairs of molecules are loaded into the cells of 
an optical potential. Whether K-K and Rb-Rb or K-Rb 
pairs are loaded together depends upon the ratio of the 
interaction strength between all pairs. Using the low- 
field resonance shown in Fig. one can clearly vary the 
interaction from strongly attractive, thus favoring for- 
mation of K-Rb pairs, to strongly repulsive, thus favor- 



3.0 



2.0 



°s 1.0 
o 



0.0 



0J 

Pi 



-1.0 



-2.0 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 

f=2 




: 


1(1,1) -KurV 






M 


i(i,-d+ (i,-i) 




f=3 










- 







B : 

/ ' " 



25 50 75 100 125 150 
B(G) 



FIG. 3: Field dependence of the purely real on, 11 and of 
the real part 7£ai_i,i_i for 41 K- 87 Rb collisions and the nom- 
inal a 3 , at, Ce- Approximate zero- field quantum numbers / 
, with f — f a + fb, are assigned to the resonances. Near the 
(2, — 2)+(l, 0) channel threshold B t (vertical line) 7?.ai_i,i-i 
changes dramatically. 



ing K-K and Rb-Rb pairs. For instance, the condition 
aKRb(u,n) = 0. 5^^(11,11)0^(11,11) required in Ref. 
in order to load a single K and a single Rb per cell can 
easily be reached. A Raman photoassociation pulse or 
a direct microwave pulse can then be used to associate 
atoms into molecules at each lattice cell. In alternative, 
a time-dependent magnetic field has also been proposed 
for near-unity conversion of atoms into molecules, see 
Ref. • This technique can of course be used for both 
the isotopic pairs investigated in this work. 
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